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£SJ . We have studied the antiferromagnetic Ising chain in a transverse magnetic field h x and uniform 

longitudinal field h z . Using the density matrix renormalization group calculation combined with a 
finite-size scaling the ground state phase diagram in (h x ,h z ) plane is determined. It is shown that 
there is an order-disordered transition line in this plane and the critical properties belong to the 
universality class of the two-dimensional Ising model. Based on the perturbation theory in h z the 
scaling behavior of the mass gap in the vicinity of the critical point (h x — 1/2, h z = 0) is established. 
It is found that the form of the transition line near the classical multicritical point [h x = 0, h z = 1) 
is linear. The connection of the considered quantum model with the quasi-one-dimensional classical 
Ising model in the magnetic field is discussed. 
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Recently, the study of the field-induced effects in low-dimensional quantum spin systems has been attracting much 
interest from theoretical and experimental points of view 0,0, For the system with an anisotropy of exchange 

interactions the magnetic properties essentially depend on the direction of the applied magnetic field. For example, 
C3 ■ the behavior of the one-dimensional antiferromagnetic XX Z model in a transverse magnetic field is drastically dif- 
ferent in comparison with the case of the longitudinal field. In particular, the transverse field induces the staggered 
magnetization in the perpendicular direction and the continuous phase transition takes places at some critical field 
This effect has been observed in quasi-one-dimensional antiferromagnet CS2C0CI4 0, where the mag- 
netic field has both the transverse and longitudinal components. Therefore, it is important to study the properties of 
the antiferromagnetic s = h XX Z model in mixed transverse and longitudinal magnetic fields 

^ : h = J2(s x n s* n+1 + sys% +l + A5*s* +1 ) - k Y,sz - /», £s* (1) 

We consider the most simple case of this model - the antiferromagnetic Ising chain given by the Hamiltonian 
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H = ^ S n S n+l - K S n - h z S n ( 2 ) 



In spite of the simple form of this Hamiltonian it cannot be solved exactly. In [lflj the ground state phase diagram 
of this model has been investigated using numerical diagonalization of finite systems and finite size scaling procedure. 
It was found 01 that the transition line between the ordered and disordered phases exists and it was assumed that 
the model belongs to the universality class of the two-dimensional Ising model. 

We note that the ferromagnetic Ising chain in the mixed fields has been studied intensively UH , , - Though 
the ferromagnetic and the antiferromagnetic Ising chains in the transverse magnetic field are equivalent, properties of 
these two models are very different at h z 7^ 0. (In fact, the model © can be transformed to the ferromagnetic chain 
but in a staggered longitudinal field.) For example, the ground state phase transition in the ferromagnetic model 
is smeared out by the longitudinal field in contrast to the antiferromagnetic model for which the phase transition 
remains at h z 7^ 0. 

In this paper we study the model (J2J) using the Density Matrix Renormalization Group (DMRG) technique 
This method allows us to consider the systems up to a few hundred sites and to determine the transition line with 
high accuracy. Using the finite-size estimation of the ground state energy and low-lying excitations we will show that 
the model (0) on the transition line is described by the conformal field theory with the central charge c = |. Besides, 
we consider the behavior of the model in the vicinity of the special points h x = i, h z — and h x = 0, h z = 1 and 
determine the form of the transition line near these points. 

The paper is organized as follows. In Sec. II we will provide a qualitative physical picture of the ground state phase 
diagram based on the classical approach. In Sec. Ill and Sec. IV the behavior of the system in the vicinity of the 
special points of the critical line will be considered. In Sec.V we will present the DMRG calculation of the critical 
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line. Sec. VI is devoted to the connection of the model with the statistical quasi-one-dimensional Ising model in 
the external magnetic field. In Summary we discuss our results. 



THE CLASSICAL APPROACH 

In order to provide a physical picture of the phase diagram of the model we use the- classical approximation, 
when spins are represented as three-dimensional vectors. The classical ground state is given by a configuration in 
which all spin vectors lie in the XZ plane with the spins on odd and even sites pointing respectively at angles ipx and 
— if2 with respect to the X axis. The classical energy of this state is 

E/N = — — sin (fx sin ip 2 ^~( cos + cos P2) ^( sm V 1 — sm (^) 

The angles ipx and ip2, minimizing this energy are solutions of the following equations 

cos (pi sin (fi2 + h z cos ipi — h x sin (pi = 

sin tpi cos if2 — h z cos <p2 — h x sin if 2 = (4) 
The solution of these equations is simple in the particular cases h z = 

<Pi = <P2 — arccos(/i x ), h x < 1 

<Pi = ¥>2 = 0, h x > 1 (5) 

and h x = 

Pi = ^2 = -, ^ < 1 

= -<P2 =2' ^ > 1 ( 6 ) 

But in general case, when h x ^ and /i z 7^ 0, the phase diagram is divided on two regions (see Fig. QJ. In the 
paramagnetic region (PM) the energy minimum is given by the configuration with ipx — —f>2 — with ip determined 
by the equation 

h z cos ip — h x sin (p — sin ip cos ip (7) 
The antiferromagnetic (AF) region is characterized by non-zero staggered magnetizations in both X and Z directions 

SZn-Sin+i = sin( )sm( ) 

oz az r<f2 - <Pl s . ,*Pl + P2s ,„s 

SL-S^+x = cos( )sm( ) (8) 

with ipi + ip2 ^= 0. 

Thus, in the AF region the magnetic fields induce a perpendicular antiferromagnetic long-range order (LRO). In 
the classical approach the value {ipx + ^2) plays the role of the LRO and it vanishes on the transition line determined 
by the equations 

h x cos tp + h z sin ip = 1 

h z cos ip — h x sin ip — sin tp cos ip (9) 

where ipx — —if2 — P>- 

The solution of Eqs.Q gives the transition line in the explicit form 



h, = 



\ll-h 2 J\l + h 2 J 3 ) (10) 



The classical phase diagram is shown on Fig. ^ The transition (critical) line separates the phase with non-zero 
AFM LRO from the phase with uniform magnetization (the paramagnetic phase). 

Of course, the classical approach does not give the correct description of the phase transition. At first, the quantum 
fluctuations shift the critical point h x = 1 to h x = | at h z = 0. Secondly, the form of the critical line is certainly 
incorrect at /i x < 1. Besides, the order parameter (tpi + ^2) vanishes on the critical line with the critical exponent 
i, which is not valid at least in the critical point h x = 5, h z = 0. Nevertheless, the fact of the generation of the 
staggered magnetizations perpendicular to the field direction is qualitatively true. 
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FIG. 1: The ground state phase diagram of the model ©. The critical line between the antiferromagnetic and paramagnetic 
states obtained from the DMRG calculation is shown by thick solid line and that in the classical approximation by thin solid 
line. 



THE CRITICAL ITF POINT 



Before the numerical determination of the critical line we consider the behavior of the model (J2J) in the vicinity of 
two special points: h x — i, h z = and h x = 0, h z = 1. This will give us the form of the critical line near these points. 

In the case h z = the model J2J reduces to the exactly solvable one-dimensional Ising model in the transverse 
field (the ITF model). This model is well studied 0] and it belongs to the universality class of the two-dimensional 
Ising model. At h x < | the ITF model is gapped and there is non-zero staggered magnetization ((— l) n S*). The 
ground state at h x < ^ is two- fold degenerated in the thermodynamic limit. At the point h x — ^ the model becomes 
gapless and the staggered magnetization vanishes with the critical exponent g. This is the transition point from the 
antiferromagnetic state to the paramagnetic gapped state. The gap vanishes linearly at h x = \. 

Now we consider the ITF model with small longitudinal magnetic field. For small magnetic field h z we rewrite the 
Hamiltonian J2J) in the form 



H = H + V 

) j ^7)' ' ,i- j " > 

■E 

After Jordan- Wigner transformation to the Fermi operators , c„ 



v = -h z ys* (ii) 



= ^p- n ^ - H^) 

j<n 

SI = \ -C+Cn (12) 

the ITF Hamiltonian Hq takes a bilinear form 

H = h x Y^ (ctc n - r^J + ~ (c+ - c„) (c+ +1 + c n+ i) (13) 
The Hamiltonian Hq commutes with the 'parity' operator 

P = exp(i7r^c+c„) (14) 

because Hq can change the number of c„ excitation by an even number only. Therefore, the space of states of Hq is 
divided on two sectors with odd (P = — 1) and even (P = 1) number of the Fermi particles c„. 



4 



The Hamiltonian Hq is diagonalized exactly 



with Fermi particles r/k and spectrum 



(15) 



which gives a gap at momentum ir 



e \ = + T + cos fc 



m 



The gap vanishes at the critical point h x = 1/2, where the spectrum becomes 



A: 



cos ■ 



(16) 



(17) 



(18) 



Below in this section we will consider the perturbation theory in V Ijllfl in the critical point h x = 1/2. 

The transition operator S z =J2 &n m V TO conserves the momentum and changes the 'parity', because as it follows 
from Eq. i|12|) it changes the number of the Fermi particles Cj (and also the number of % particles) by odd number. 
Therefore, the non-zero matrix elements with the transition operator S z have the states with equal momentum but 
different parity. The last fact means that the perturbation theory in V contains only even orders. 

The second order correction to the ground state energy is: 



SE, 



(2) 



\S z \s) (s\S z \0) 
En — E s 



(19) 



The ground state |0) has the momentum q = and zero number of the rjk particles (P = 1). Therefore, the non-zero 
contribution to the sum in (|19|) is given by the intermediate states \s) with zero momentum and the parity P = — 1. 
As follows from Ea. (|18|l . all states \s) with momentum q = and odd number of the rjk particles (P = —1) have 
'high' energies E s — Eq = e s 1. On the contrary, among the states with q = and even number of the r/k particles 
(P = 1) there are many states like Vx-k f l-Tr+k 1^) with small k having small excitation energies e s ~ k and they can 
lead to infrared divergencies. 

Hereinafter we consider large but finite systems of length N. We shall study the dependence of dominant contribu- 
tions to the perturbation theory on N, omitting numerical factors. Using the fact that e s > 1 one can rewrite Ea. i|19|) 
as 



5E, 



(2) 



-IE 



\s z 



(a\S'\Q) 



(sr\o) 



(20) 



In higher orders of the perturbation series for the ground state energy each second intermediate state \s) has odd 
number of the Fermi particles (P = —1), and, therefore, high energy e s > 1. For example, let us consider the fourth 
order correction to the ground state energy 



SE. 



(4) 



s,s f ,s" 



-SE^hlY. 



\S z \s) {s\S z \s') (s'\S z \s") (s"\S z \0) 
{E - E.) {E - E s ,) (E - E.„) 

S z \s) {s\S z \Q) 



(E - E s y 



(21) 



All intermediate states s,s',s" have momentum q — 0. The states s' have P = 1 and some of them have small 
excitation energies e s i ~ l/N, while the states s and s" have P = — 1 and high energies e s > 1. Therefore, one can 
sum over the intermediate states s and s" , which reduces Ea. (|21|l to 



SE, 



(4) 



(0\(S z ) 2 \s>) (s'\(S z ) 2 \0) 

En — E s i 



(22) 
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Now we note, that this expression looks like the second order correction and Eq.J3U|) looks like the first order 
correction to the ground state energy with a perturbation — h z (S z ) . In a similar way one can sum over all intermediate 
states with q = and P = — 1 in all orders of perturbation series. As a result, we arrive at the perturbation theory 
with the effective perturbation 

V 1 = -h 2 z (S*) 2 = -h 2 z J2 S nS Z m (23) 

n, m 

We note, that the perturbation theory with the perturbation V\ coincides with the original perturbation theory 
(|II|I in a sense, that both perturbation series have the same order of divergencies (or power of N) at each order in 
h z . But numerical factors at each order in h z can be different. 

The perturbation V\ commutes with the parity operator Ijl4(l and conserves the momentum. Therefore, the pertur- 
bation series in V\ contains the intermediate states with q = and P = 1 only, and some of these states have small 
excitation energies e s ~ 1/N. These states give the most divergent contribution to the perturbation series in V\ and 
further we shall take into account these states only. 

Now we need to estimate the matrix elements of the operator (S z ) 2 . The behavior of the correlation function 
jffijj iS 1 ^) in the ground state and in the low-lying states with excitation energies e s ~ 1/N on large distances is known 

{S'nS'J ~ ^ ' |1/4 (24) 
\n — m\ 

and, therefore, due to oscillation of the correlator (S^S^) the sum over n and m can be estimated as 

E < 5 » 5 m> = 1[ + 2N J2 ( S ^n) - «-07465(l)iV (25) 

n,m n>l 

where the constant was found from extrapolation of the exact results for finite chains with N = 6, . . . 14. 

The non-diagonal matrix elements of the operator (S z ) 2 with two different low- lying states was also calculated 
numerically. It was found that the only non-zero matrix elements (in the thermodynamic limit) are given by the 
states s,s' differing by two rjk particles, like \s') — rf^Ti^ k \s). All such pairs of states give the same value for the 
matrix element 

(s| (S**) 2 = 0.3108(1) (26) 

while all other matrix elements exponentially drop with N. 

Now we estimate all terms of perturbation series in h z . According to Ea. l|25|l the second order correction l|20|l is 
proportional to h 2 N. The factor at h 2 N was found numerically by exact diagonalization of finite systems and the 
following calculation of the sum (|19f) . The data for N = 6, ... 14 are well extrapolated and give 

5E {2) = -0.07060(5)/i 2 iV (27) 

According to Ea.l|27|l the zero-field susceptibility is Xz = 0.1412(1). This result is in a perfect agreement with the 
value Xz — 0.14118 . . . obtained analytically in |l8j |. 

Using Ea. (|26|l we find that the non-zero contribution to the fourth order correction to the ground state energy 
(1221) is given by the states |s) = r]^rj_ k |0) with the excitation energy s s = 2 cos 4. Thus, because of small excitation 
energies at k ~ it in denominator in ll'L'l) the fourth order correction to the ground state energy turns out to be 
proportional to TV (we omit here logarithmic corrections) 



J2r~ h * N (28) 



The ra-th order in V\ is proportional to h 2m . The denominator of the m-th order contains (m — 1) small excitation 
energies ~ 1/N and all matrix elements in the numerator are of the order of unity. Therefore, the m-th order in V\ 
diverges as h 2m N m . But this is not valid for odd orders in V\. The analysis shows that all odd orders in V\ diverge 
as h 2m N m ~ 2 . For example, the third order correction in V\ (the sixth order in h z ) for the ground state energy does 
not diverge: 

SE^ - h\N (29) 
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So, the odd orders in V± give the next order corrections to the ground state energy and we omit it below. 
Summarize all above, we arrive at the perturbation series in the form 

oo 

5Eo = -a h 2 z N-h' l z NY J b n (h 2 z N) n (30) 

with ao — 0.07059 . . . (Ea.l|27|l) and unknown constants b n . One can see that the first divergence appears only in the 
eight-th order in h z , and it is very difficult to observe it numerically. The exact numerical calculations up to N = 14 
of the second, the fourth and the sixth order corrections in h z confirm the form of the series 1)30(1 . 

The sum in Ea. (|30|l forms the scaling function fo (x) of the scaling parameter x — h 2 N. Thus, the ground state 
energy takes the form 

5E = -a Q h 2 z N - h 4 z Nf a (x) (31) 

Since the ground state energy is proportional to N, the scaling function fo (x) has finite thermodynamic limit at 
x — > oo. Thus, the leading term of the perturbation theory for the ground state energy is given by the second order 
and the divergent part of the perturbation theory gives the correction ~ h\. 

The perturbation series for low- lying states has the same form as in Ea. (|31|l . but each low- lying state has its own 
scaling function. Therefore, for the first excited state (with momentum q = it) one has 

SE„ = -a^h 2 z N - h z NU (x) (32) 



So, the mass gap m = SE n — SE appears as: 

m — a m h 2 + g (x) h 2 (33) 

where a m = (ao — a„)N and the scaling function g (x) = x (/ (x) — /„. (x)). Since the gap is finite, the scaling function 
g (x) in the thermodynamic limit at x — > oo must tend to some finite limit 

m = a m hl + g (oo) h\ (34) 

From the last equation we see that the gap is proportional to h 2 , but the factor at h 2 is given not only by the 
second order correction a m but by all collected divergent orders of the perturbation series. The numerical estimation 
of the second order correction to the gap gives 

a m = 0.1875(2) (35) 

In order to find the factor at h 2 in mass gap we performed DMRG calculations (for details see Sec.V) of the model 
(|ll|l for h z = 0, . . . 0.3 and N = 20, ...300. The dependence of the gap on the scaling parameter x = h 2 N is shown on 
Fig. [3 One can see that the points with different N and h z lie perfectly on one curve. It definitely manifests that 
the scaling parameter is h 2 N . The DMRG calculation data give the mass gap 

m = 0.37(1)^ (36) 

From this equation follows that the second order correction to the mass gap (|35|) gives approximately one half of a 
total gap, and another half is collected by all other divergent orders contained in the scaling function g (x). 

We note, that in contrast to Ea. (|36|l for the ferromagnetic ITF model the mass gap at the critical point h x = \ is 

proportional to h^ 15 (l2^.[l3j. 



Mapping to the effective model 

Let us return to the estimation of the diagonal element of the operator V\ (|25|) . Using the fact of the oscillation of 
the correlator l|24|l in the ground state and in the low-lying states we can rewrite the sum in (|25f) approximately as a 
one half of the first term: 

(^) 2 ) = J + 2 E ( S nSrn) * J + £ (SfAI+l) (37) 
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FIG. 2: Scaled mass gap S f° r h x = i with various values of N and h z as a function of the reciprocal scaled parameter —j-. 



The last equation suggests that the perturbation V\ can be reduced to the operator V% 

Vx^V 2 = -ah 2 z N - bhl Y, S z n S z n+1 



(38) 



with some constants a, b, which we will define later. In order to verify this assumption we compared the matrix elements 
of the operator V2 with those of V\ . We have found that the dependence of the matrix elements of the operator V2 
on N is in a full accord with Eas. (|25[) and i|26|) . That is, the diagonal matrix elements of V2 are proportional to N 
and the only non-zero (non-diagonal) matrix elements have the pairs of states differing by two rjk particles 



(s| S*S* +1 \s) — - 

Thus, we arrive at the effective Hamiltonian 

H cfi = -ahlN + (1 - bhl) E S n S n+i -\T. S n 



(39) 



(40) 



Again, the original model (f 1 1|) is equivalent to the effective model (1401) in a sense, that the perturbation series for 
both models have the same order of divergencies (or power of TV) at each order of h z . But numerical factors at each 
order of h z can be different. 

The advantage of this mapping lies in the fact that the effective model l|40() is of the ITF type and, therefore, is 
exactly solvable one. The spectrum of H e g is 



2 1 + J , J u 
e k = — h — cos k 

4 2 



(41) 



with J = 1 — bh 2 z . So, according to the spectrum of the effective model (|41|l the magnetic field h z produces the gap 
in the spectrum 



2 z 



(42) 
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As was noted above, apart from numerical factors the effective model (|40|l has the same form of perturbation series 
in h z as the original model (|f lfl . Therefore, we can check the form of the perturbation theory (|3(JI) by studying the 
perturbation series in h z of the exactly solvable model (|40|) . 

For large but finite systems the ground state energy (with q = 0) and the first excited state (with q = ir) of the 
Hamiltonian l|4U|l are known exactly (lij: 



Eo = -^-ig^ + 4(l- 9 )si^(|-^J 



7T \ 



K = -a^JV+f-ix! X 2 +4(1 - 5 ) sin 2 f^) (43) 



4 4^V ^ \N 

j=i V v - 

where g = 

The perturbation series for the ground state energy and the gap can be easily found by formal expansion of these 
expressions in small parameter g. This results in the scaling form of perturbation series for the mass gap: 



m= 9 j+g 2 NY;dn{gN) 2n (44) 



■" • -^VV <..\-> 2 " 
ri=0 

where constants d n determine the scaling function G(x) of a scaling parameter x = gN 



m= 9 -+gG{x) (45) 

The scaling function G(x) was found in |19|. In the thermodynamic limit x — > oo the function G(x) — > 1/4, resulting 
in the gap (|42fl . One can see that the scaling parameter and the form of the perturbation series 144|) coincides with 
Eq.(|2Sl- 

Now we define the constants a and b in Eg. (|40fl . so that the second order corrections to the ground state energy 
and the gap of the effective model (|4T1)l numerically coincide with Eas. f?7jl and This leads to the values 

a = 0.190(1) 

b = 0.750(1) (46) 
The mass gap of the effective model Ij40ll with constants a and b defined above is 

m = 0.375(l)/i 2 (47) 

Surprisingly, the value of the mass gap of the effective model turns out to be very close to one found numerically 
for the original model (Eq.J2U). It means, that the gap scaling function G{x) of the effective model (|40|l has the 
same thermodynamic limit as that for the model (|ll(l . Moreover, we have calculated numerically the fourth order 
corrections in h z to the ground state energy and the gap for the models (|11J) and (|40(l and found perfect numerical 
agreement (the second order corrections in h z for both models coincide by the definition (|46Jl ). Therefore, we expect 
that the mapping of the model (|ll|l at h x = 1 /2 to the model l|4U|l with constants a and b defined in (|46|l is exact for 
low- lying excitations in the thermodynamic limit. 

The mapping of the model (jll|) to the effective ITF model can be extended to the case h x ^ 1/2. For this case the 
effective Hamiltonian becomes 

#cff = ~a(h x )h 2 z N + (1 - b(h x )h 2 z ) S z n S z n+1 ~KJ2 S% (48) 

where now the constants a, b are the functions of h x . The functions a{h x ) and b(h x ) are defined in such a way that 
the second order corrections to the ground state energy and the gap for the effective model (|40f) coincide with those 
for the original model (fill) . 

Using the effective Hamiltonian (|48|l one can calculate the susceptibility Xz(h x ). In the vicinity of the critical point 
h x = 1/2 it has a logarithmic singularity 

Xz(h x ) = Xzc ~ ^ \h x - log U x - i J (49) 



with Xzc = 0.1412 from Eq.^ and b = 0.75 (Eq.gnj). 
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THE MULTICRITICAL POINT 



Another exactly solvable limit of the model J5J) is the case h x = 0. The model at h x — is the classical one. At 
h z < 1 its ground state is the antiferromagnet. At h z = 1 the first order phase transition to the ferromagnetic ground 
state occurs. This is so-called multicritical point. The Hamiltonian at h~ = 1 has a form 



H = H Q + V 

H o = -j + ( 5 « ~ 2) ( 5 ' n + l 2 
V = -^(S+ + S-) (50) 



where S± = E#n- 

The ground state of Hq is macroscopic degenerate: all spin configurations, excluding those with two neighbor spins 
pointing down, have the same energy — x- The number of these states is ( 1 + 2 v/ ^ ) JV |2(ij . The transverse field h x lifts 
the degeneracy. The exact calculation of the first order correction in h x for N ^> 1 is rather complicated because 
it involves the exponentially large number of degenerate states. We carried out the approximate calculation of the 
perturbation theory in V within degenerate manifold using a simple variational function in the form 

N/2 
m=0 

where \l/ m are the sum (with equal weights) of all admissible states with m spins down: 

Vm=w~ i (S-p) m * F (52) 

Here P is a projector excluding states with two neighbor down spins, ^Sfp is the ferromagnetic state with all spins 
up and the normalization factors w m are 

(N-m-l)\N 

Wm = m\(N-2m)\ (53) 
Matrix elements of V with respect to \P TO are 

{^mV^m') = Y(«m^m',m-1 + %+l^m',m+l) (54) 

where 



N — m 

Coefficients c m in (|51|l obey the equations 



(N -2m + 2)(N -2m+l) 

(55) 



(E + -j)Cm = y(Sm+lC m +l + 8 m Cm-l) (56) 



The quantity s m has a sharp maximum for JV > 1 at m = 3 N and s mo is 



s mo = y(\/5-l)\A/5-2 (57) 

The ground state energy in the thermodynamic limit is defined as a lowest eigenvalue of Eas. (l56|) . which is 

TV AT 

So = - j ~ h x s mo =-—- 0.30028^/1* (58) 
In the frame of variational function approach H51JI the ground state magnetizations (S^) and (S*) are 

= = 0.300028 (59) 



10 



The function i|51|) has a momentum q — 0. To calculate the spectrum in the first order in h x it is necessary to 
choose the variational function of the type (|51|l with momentum q. We omit here rather cumbersome calculations, 
which shows that the mass gap corresponds to q — n and 

m = h x (60) 

The numerical diagonalization of finite cyclic systems shows very rapid exponential convergence to the thermody- 
namic limit. The extrapolated values for the ground state energy and the mass gap (at q = n) for h x <C 1 are 

§ = -\ -0.3017(1)/** 
m = 0.4841(1)/^ (61) 

Comparison of these results with Eg. (|58|l shows that the variational ground state energy differs from the 'exact' 
one within 0.4%. As follows from Eq. l|61|l the mass gap in the multicritical point opens linearly with h x , which is 
correctly described by the variational approach (|60(l . 



THE CRITICAL LINE 



In general, the critical line at h x ^ and h z ^ can not be found exactly. To obtain it we used the DMRG technique 
(l4|. We have performed DMRG calculations using both the infinite-size and the finite-size DMRG algorithms. We 
calculated the ground state energy E (N) and two lowest excitations mi(JV) = Ei(N) — E (N) and m,2(N) = 
E 2 (N)-E (N). 

In order to check the accuracy of the DMRG method we compared the obtained results with the exact ones for 
the ITF model. We used the infinite-size algorithm and open boundary conditions. The dependence of the results on 
the number of retained states s in the DMRG computation and on a number Nrg of DMRG steps (N = 2Nrg + 2) 
has been investigated. We have found that the calculation with s — 25 gives satisfactory accuracy up to N = 300. 
For example, relative errors in the ground state energy and in the mass gap at h x = ^ are — 10~ 9 (10~ 7 ) and 
= 10~ 7 (10~ 5 ) for N = 100(300). The accuracy becomes better when the value \h x — ^| is increased. 

The critical field h xc (h z ) at a fixed value h z (0 < h z < 1) is determined by vanishing of the gaps mi and mi. Below 
the critical field the mass gap mi (TV) — > exponentially with N. This behavior confirms the fact that the ground 
state is doubly degenerate in the thermodynamic limit at h x < h xc (h z ). The true mass gap in this region is defined 
by the value of mi- The typical behavior of the gaps miand mi extrapolated to the thermodynamic limit is shown 
on Fig. |3| We note that at h x > h xc (h z ) the gaps miand m 2 coincide at N — > oo. 

The critical line obtained from the DMRG results extrapolated to the thermodynamic limit is shown on Fig.^ In the 
plane (h x , h z ) it connects two limiting critical points studied in previous sections: the critical point (h x = g, h z = 0) 
of the ITF model and the multicritical point (h x = 0, h z = 1). Now we discuss the properties of the model J5J) near 
the critical line. 

The ITF model in the critical point is described by a conformal field theory with a central charge c = i. We 
expect that the ITF model is generic for the model (0) on the whole critical line (except the multicritical point 
h x — 0, h z = 1). To verify this suggestion we estimated the value of c on the critical line. For this we used well 
known fact |2ll | that for conformal invariant model with periodic boundary conditions the central charge appears at 
1/7V correction to the ground state energy 

7TCV 

E = e 00 N~— (62) 
6iV 

where eoo is the ground state energy per site in the thermodynamic limit and v is the sound velocity. 

At first, we calculated the sound velocity v as (similar to the spectrum of the ITF model l|18J) the sound velocity 
on the whole critical line is determined at k = n) 



N 

v(N) = - 

Z7T 



2-7T 

E(n--)-E(n) 



(63) 



We carried out these calculations using the numerical diagonalization of the Hamiltonian with the periodic 
boundary condition for N < 14. The size extrapolation is carried out by formula v(N) = v + aN~ 2 . The example of 
the extrapolation procedure is shown on Fig. The dependence of v on h z along the critical line is shown on Fig. [S] 
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hx 

FIG. 3: The mass gaps mi (triangles) and 7712 (circles) as functions of the transverse field h x at h z = 0.9 from the DMRG 
calculations extrapolated to the thermodynamic limit. 




0.470 1 1 1 1 1 1 1 

0.000 0.005 0.010 0.015 0.020 0.025 0.030 

1/N 2 

FIG. 4: The extrapolation procedure of finite size dependence of the sound velocity in the critical point h x — 0.4962, h z — 0.1. 
(v = 0.4960). 

After that, the central charge c has been calculated using Ea. (|62(l . The size extrapolation shows that the central 
charge is c = 0.500(1) for all calculated critical points. 

The critical exponents v of the ground state correlation function (OiOr) ~ R~ U (R ^> 1) is defined by the scaling 
dimension X of the operator O and v = 2X [22| • The scaling dimension is related to the finite-size correction of the 
lowest excited eigenstate which can be reached from the ground state by the operator O. The finite-size correction to 
this excitation energy is 

2ttv 1 

AEi = E Z -E = —X t + o(-) (64) 
For the ITF model at h x = \ the sound velocity is v — \ . The ground state has a momentum q = and the lowest 
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FIG. 5: The dependence of the sound velocity on h z along the critical line. 



excitation energy in the sector with momentum q = is 



and in the sector with q = n is 



E(q = 0)- E = 2 sin — 



E(q = tt) - E Q = - tan ^ 



(65) 



(66) 



Thus, the scaling dimensions are Xq = 1 and X T = |. The corresponding associated operators are 

n 

o q ^ = £(-i) n 5* 

n 

which is in accord with the well known results [l7| for the asymptotic of the correlation functions 



(67) 



W0 J -R/ WO/ 

(S$S} 



R- 2 
R 1 /* 



(68) 



Using Eq. I|64|) we numerically checked that the scaling dimensions related to the finite-size corrections of the lowest 
excited eigenstates on the critical line remain as in ITF model X Q = 1 and X w = §. Therefore, we conclude that 
the model (J2J on the critical line belongs to the universality class of the ITF model. This means that in accord with 
the prediction of the classical approach to the left of the critical line the staggered magnetizations ((— l) n 5*) and 
((— l) n S%) exist. But in contrast to the classical approach they vanish on the critical line with the critical exponent 
g. The mass gap is closed on the critical line and the critical exponent for the gap is equal to unity, i.e. the scaling 
behavior of the gap near the critical line is linear. For the particular case h z = 0.5, this is illustrated on Fig.0 where 
the scaling plot for the scaled mass Nm\ with N{h x — h xc ) is shown. 

The behavior of the critical line near the ITF point h x = 1/2, h z = can be found from the following consideration. 
As it was established above in the vicinity of the critical line the gap is proportional to deviation from the line. This 
is valid for any direction of deviation except the direction at a tangent to the critical line. In the vicinity of the ITF 
critical point for fixed according to Ea. (|17|l the gap is m = \h x — h xc (h z )\. On the other hand, the gap is given 

by Ea . IJ36JI ■ Combining these two expressions for the gap on the line h x = 1/2, we obtain equation for the critical line 
in the vicinity of the point h x — 1/2, h z — as 



0.37/1? 



(69) 
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FIG. 6: Scaled mass gap Nm\ for h z — 0.5 with various values of N plotted as a function of the scaled field Nh x . 



This expression for the critical line is in a very good agreement with our numerical estimations up to h z ~ 0.5. 

In the vicinity of the multicritical point the form of the critical line can be found in a similar way. At h z — 1 and 
h x <C 1 the mass gap is proportional to h x (see Ea. ()f>ll) V On the other hand, for fixed h x <C 1 this gap is proportional 
to a deviation from the critical line m ~ (h z — h zc ). Therefore, the critical line near the multicritical point behaves 
as h xc = A(l — h z ). The numerical coefficient A has been found from the DMRG calculation of the critical line in the 
vicinity of the multicritical point. As a result, the form of the critical line at h x <C 1 is 

h xc (h z ) = 1.50 (1 -h z ) + ((1 - ^) 2 ) (70) 
This form differs from that given in [lo|. where h xc (h z ) ~ (1 — h z ) ' 75 . 



RELATION WITH THE TWO-DIMENSIONAL ISING MODEL 



It is well-known that the critical properties of the two-dimensional Ising model are equivalent to those of the one- 
dimensional ITF model [5^ . Using the formalism of Ref.jl^l one can show that the quantum Hamiltonian (J2J is 
related to the transfer matrix of the strong anisotropic (quasi-one-dimensional) Ising model in the uniform magnetic 
field. This model describes ferromagnetic Ising chains weakly antiferromagnetically coupled with each other. The 
Hamiltonian of the model is 

H = —J\ CTn,mCn+l,m + J2 J]] Vn,m&n,m+1 — h <Jn,m (71) 

where cr n>m = ±1, J\ ^> J2 > 0. 

The relation between the model © and the model (|71|l are given by 

2h x f3J 2 = e~ 2f3Jl < 1 

h = 2h z J 2 (72) 

with 13 = 1/kT. 

The equivalence of these two models means that the model 11710 undergoes the phase transition from the ordered (at 
T < T c ) to the disordered (at T > T c ) phase. The critical line T c {h) of the quasi-one dimensional model are related 
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to the critical line h xc (h z ) of the model By use of mapping (|72() we obtain that the critical temperature T c (h) is 

2 7 

kT c (h) = — ^ (73) 

l0g \J 2 h x (h/2J 2 ) J 



In particular, at h — ► 



and T c — > when h — > 2 J? as 



fcT c = (74) 



fcT c = J J \ N (75) 



log 



The order parameter i (cr n ,m + 0Vi,m+i) vanishes at T — > T c with the critical exponent 1/8. We note that the phase 
transition does not occur if J2 < as well as in the ferromagnetic version of the model J5J. 

The free energy of the model l|71|) is related to the ground state energy of the Hamiltonian by 

F = AJ 2 E (h x ,h z ) (76) 

According to Eas. l|49l) . (|72|l and (|76[) the zero-field susceptibility in the vicinity of the critical temperature T c (0) 
can be obtained using exactly solvable effective model l|4UI) . At \T — T c (0)| <C T C (Q) the susceptibility is 

0.1412 -0.24 * ™1„ 6 (™>) ,77, 



fcr c (o) t c (o) °V r c (o) 

The last equation has a form coinciding with the results obtained for the two-dimensional antiferromagnetic Ising 
model by the series expansion method |24| . 



SUMMARY 



We have studied the antiferromagnetic Ising chain in the mixed transverse and longitudinal magnetic field. It 
was shown that the quantum phase transition existing in the ITF model remains in the presence of the uniform 
longitudinal field. Using the DMRG simulations we have found the critical line in the (h x , h z ) plane where the mass 
gap is closed and the staggered magnetizations along the X and Z axes vanish. It is found numerically that the model 
on the critical line is described by the conformal field theory with the central charge c = 1/2, i.e. it belongs to the 
universality class of the ITF model. 

The scaling behavior in the vicinity of the ITF critical point is studied in detail. It is shown that the mass gap 
is proportional to h 2 z and the contributions to it are given not only by the second order correction but also by all 
other divergent orders of the perturbation series in h z . Besides, the analysis of the perturbation theory shows that 
the considered model at h z <g; 1 can be mapped to the effective ITF model with renormalized parameters depending 
on h z . In a framework of the effective ITF model the behavior of the susceptibility at h x ~ 1/2 is determined. 

The behavior of the model in the vicinity of the multicritical point is investigated. Using both the variational 
approach and the numerical diagonalization results we have found that the mass gap is proportional to h x . Close to 
the multicritical point the form of the critical line is linear. 

Of course, the considered model is the simplest case of the XX Z model It is interesting to extend the present 
analysis to this model as well as to study effects of interchain interactions in quasi-one-dimensional generalization of 
the model ©• 

The mapping of the quantum model (0) to the strongly anisotropic statistical two-dimensional Ising model in the 
uniform magnetic field was considered. The behavior of the susceptibility of this model near the critical temperature is 
found. We expect also, that the main features of the considered model are valid for the statistical two-dimensional 
antiferromagnetic Ising model in the uniform magnetic field h. That is, the applied magnetic field does not smear the 
phase transition existing in the two-dimensional Ising model, which is generic case for the whole transition line T c (h) 
in the plane (T, h). In particular, we believe that the analysis of the perturbation series in h <C 1 in the same manner 
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as was done in Sec. Ill gives the scaling behavior for the correlation length and the form of the critical line near the 
point T c (0) is T c (0) - T c (h) ~ h 2 . 

This work is supported under RFBR Grant No 03-03-32141 and ISTC No 2207. 
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